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NONLINEAR  MODE  COUPLING  THEORY  OF  THE 
LOWER-HYBRID-DRIFT  INSTABILITY 

1*  Introduction 

The  lower-hybr id-d rift  instability  has  been  of  considerable  Interest 
to  plasna  physicists  for  more  than  a  decade  as  a  driving  mechanism  for  the 
anomalous  transport  of  particles,  momentum  and  energy  in  both  laboratory 
and  space  plasmas.  In  laboratory  applications,  this  instability  was 
initially  studied  to  explain  anomalous  sheath  broadening  in  theta  pinch 
implosions,  but  has  recently  been  applied  to  other  laboratory  confinement 
devices  such  as  toroidal  reversed  field  pinches^  and  compact  torii.^  In 
space  plasma,  it  has  been  proposed  as  a  mechanism  to  provide  anomalous 
resistivity  for  reconnection  events  in  the  magnetosphere.4  Experimentally, 
the  lower-hybrid-drift  instability  has  been  observed  in  a  recent  theta 
pinch  experiment,^  and  satellite  data  Indicate  that  it  exists  in  the 
earth's  magnetotail  and  magnetopause. * 

The  linear  theory  of  the  lower-hybrid-drift  instability  is  well 

7_Q 

understood.  7  Ths  mode  can  be  excited  In  inhomogeneous  plasmas 

1/4 

wham  Pj/Lg  >  (m#/m^)  where  p^  is  the  mean  ion  Larmor  radius  and  1^  is 

the  acale  length  of  the  density  gradient.  At  maximum  growth  it  is 

characterised  by  «  <  «...  y  <  w  ,  k*B  -  0,  and  kp  >•  1  where 

r  ^  id  v  r  v  v  es  «  #  m 

w  ■  er  +  iy,  eo^  is  the  lower  hybrid  frequency  and  pft#  ■  p^m^/Zm^) 

Since  the  instability  is  high  frequency  (u  »  ft^)  and  short  wavelength 
(kp^  »  1)  the  ions  are  treated  as  unnagnetised  while  the  electrons  are 
treated  as  strongly  magnetised.  In  the  weak  drift  regime  <  1  (where 

V££  is  the  ion  diamagnetic  drift  velocity  and  v^  is  the  ion  thermal 
velocity),  o)  <  kV^  and  the  mode  is  driven  unstable  by  the  resonant 
interaction  of  the  drift  wave  with  the  ions.  In  finite  g  plasmas, 
electrons  also  resonate  with  the  wave  via  their  VB  drift.  The  resultant 
electron  damping  is  stabilising  for  kp  >  1.  In  the  high  drift 
regime  (V^/v^  >  1),  the  mode  is  driven  unstable  via  the  interaction  of  a 
positive  energy  lower  hybrid  wave  and  a  negative  energy  drift  wave.  In 
this  paper  we  focus  on  the  weak  drift  regime. 

The  nonlinear  development  of  the  lower-hybrid-drift  instability  is 
complex  and  not  as  well  understood  as  the  linear  theory.  A  variety  of 
nonlinear  saturation  mechanisms  have  been  proposed  to  date.  Ion  trapping^ 
and  stochastic  electron  heating^ *  have  been  proposed  to  quench  the  growth 
of  the  instability  in  particle  simulations.  However,  ion  trapping  is  not  a 
Haase eripl  approved  September  7,  1988.  « 


12 

viable  saturation  mechanism  when  a  broad  2D  spectrum  of  waves  is  excited, 

and  it  has  not  been  observed  in  recent  simulations  using  realistic  mass 

13 

ratios  where  such  spectra  develop*  Stochastic  electron  heating  does  not 

onset  until  a  rather  large  threshold,  n/n  >  0.25,  is  exceeded  where  n  is 

the  density  fluctuation  of  the  wave.  A  wave  energy  bound  based  on  the 

available  free-energy  in  the  relative  electron-ion  flow  (current)  has  also 

14 

been  invoked  to  calculate  a  maximum  wave  amplitude.  In  a  finite  8 

plasma,  however,  the  current  and  self-consistent  magnetic  field  are  coupled 

and  the  magnetic  field  energy  can  also  be  tapped. ^  The  free  energy  then 

becomes  extremely  large  and  does  not  realistically  act  as  a  bound. 

Finally,  electron  resonance  broadening  has  been  proposed  as  a  saturation 

mechanism. It  has  been  recently  shown  in  more  refined  calculations  that 

"resonance  broadening"  does  not  cause  a  net  dissipation  of  energy  but 

merely  leads  to  an  exchange  of  energy  between  modes  in  k  space  unless  wave- 

17  18  ^ 

particle  resonances  are  included.  ’  This  mechanism  can  therefore  only 

10 

stabilize  the  wave  spectrum  in  finite  0  plasmas. 

In  this  paper  we  present  a  nonlinear  mode  coupling  theory  of  the 

lower-hybrid-drift  instability.  Our  preliminary  calculations  have  been 

20 

presented  previously.  The  basic  result  is  that  the  instability  can 
saturate  by  transferring  energy  from  growing,  long-wavelength  modes  to 
damped,  short-wavelength  modes.  This  saturation  mechanism  is  consistent 
with  recent  computer  simulations  and  with  experimental  observations.  In 
this  paper,  we  extend,  the  earlier  theory  to  self-consistently  incorporate 
finite  B  effects,  and  to  include  electron  collisions  in  order  to  compare 
our  results  with  experimental  observations. 

The  organization  of  the  paper  is  as  follows.  In  the  next  section  we 
derive  the  nonlinear  wave  equation  and  the  wave  energy  transfer  with  the 
two  basic  nonlinearities:  the  nonlinear  electron  E  x  B  and  polarization 
drifts.  In  Section  III  we  discuss  the  numerical  methods  used  to  study  the 
mode  coupling  process.  In  Section  IV  we  present  results  based  upon  VB 
electron  damping  of  the  short  wavelength  modes,  while  in  Section  V  we 
present  results  based  upon  collisional  damping  of  the  modes.  In  Section  VI 
we  summarize  our  theory,  compare  our  theoretical  results  with  experimental 
observations,  and  discuss  some  applications  of  our  results  to  space 
plasmas . 


II.  Derivation  of  the  Mode  Coupling  Equation 


Ve  consider  a  slab  equilibrium  of  cold  electron  and  warm  ions  with  a 

A 

density  profile  n^(x)  supported  by  a  magnetic  field  B  -  B0(x)ez  as  shown  in 
Fig.  1.  Equilibrium  pressure  balance  requires 

fj  [PQ(x)  +  (x)/8ir  ]  -  0.  (1) 

2 

The  ratio  of  thermal  to  magnetic  pressure,  B  “  SimT^/B^,  is  assumed  to  be 
of  order  unity  so  that  electromagnetic  corrections  to  the  lower-hybrid- 
drift  instability  must  be  retained.  The  equilibrium  ion  velocity  distri¬ 
bution  is  taken  to  be  a  Maxwellian  with  an  average  drift  velocity  ■ 

Vdi  ey  where  Vdi  »  (v2/2nt  )3  in  n/3x,  v2  -  2^/n^  and  -  eB^/m^c. 

Because  of  the  flute  nature  of  the  instability  (k*B  ■  0),  we  self- 
conslstently  limit  the  spatial  variations  to  the  xy  plane.  In  this  limit 
the  electric  and  magnetic  fields  can  be  represented  by  the  scalar  potential 

O 

4>  and  the  vector  potentials  A^  and  Ay  as 

B  »  B  e  ■  (3A  /«>x  -  3A  /a  )  e  (2) 

-  a  *  '  y  x  y-'x 

E  -  -  76  -  c”1 3  A/at  (3) 

^  T  V 

with  V«A  -  0.  Even  when  B  ^  1>  the  induction  component  of  E  is  small  so 
that  it  can  generally  be  discarded  except  when  the  V$  component  of  E  does 
not  contribute. 

We  are  Interested  in  time  scales  of  order  9 /at  >>  so  that  the  ion 
response  to  both  the  equilibrium  and  perturbed  magnetic  fields  can  be 
neglected.  Furthermore,  in  the  weak  drift  regime  the  ions  respond  to  the 

g 

perturbed  potential  A  adiabatically  to  lowest  order  since  3/at  «  v^lvl. 
Thus, 

ni  “  +  X  Vi1|v’1|(|t  +  Vdi  1y)]exp('  (4) 

where  the  term  proportional  to  Ti  in  (4)  is  a  small  correction  describing 
the  resonant  ion  interaction. 


In  contrast  to  the  ions,  the  electrons  are  strongly  magnetised  since 
3/3t  «  Q#.  The  electron  motion  is  simply  given  by  the  E  x  B  and 
polarisation  drifts, 


.  E 

C  C  Q  v 

l.  ~Tlxtz  -  a71i  5" 

Z  6  Z 


where  d/dt  •  3/3t  +  V_*V  and  V.  •  cE  x  e  /B  .  The  electron  density  can  be 

'C  1  Xi  v  Z  Z 

calculated  from  the  continuity  equation, 

d  in  n  /dt  +  V*V  -  0.  (6) 

©  ^*© 

The  electron  compression  V»V^  is  obtained  from  Eq.  (5).  To  lowest  order, 

VVXe  ‘  bV  If  ~  It  tn  V  (7) 

z  © 

The  first  term  in  Eq.  (7)  arises  from  the  polarisation  drift  and  the  second 
from  the  compression  of  the  flux  tube  in  a  finite  B  plasma.  The  continuity 
equation  can  now  be  written  as 

If  “(V.J+rirlt  72*-°  <8) 

s  e 

so  that  the  rate  of  change  of  the  number  of  electrons  in  a  flux  tube  (n  /B  ) 

e  s 

is  given  by  the  polarisation  drift*  The  time  variation  of  Bz  in  Eq.  (8  ) 
is  calculated  from  Ampere's  Lmr 

V*Ax  -  -  <4w/c)Jx  •  -  (4n»#/Bl)3#/3y,  (9) 

where  J_  is  dominated  by  the  E  x  B  drift  of  the  electrons.  Using  V*A  -  0, 

X  v  w  v 

we  find  that  V2A_  »  -  3B  /3y  so  that 
z  z 


h  til  *  bo**)  ■  °> 


which  can  be  integrated  to  evaluate 


If  “  ■*  - 1  lr  (ff ) 


TOW 


in  Eq.  (7).  The  convective  derivative  o£  Bz  can  similarly  be  evaluated 
using  Eq.  (10)  in  conjunction  with  the  equilibrium  pressure  balance 
relation. 


«  6  c  8  tan  H 

V7  in  Bz  -  Y  B"~  ay  * 

Z 


(12) 


Finally,  invoking  charge  neutrality  (ne  -  )  and  combining  (8),  (11)  and 

(12),  we  obtain  the  nonlinear  equation 


(i  -  ?)iT  +  ;y  -  y0iv“ii(iTT  +  iTy)  + 

+  (vi  X  +  y0yj  X  e2.v|v_1|(iT  +  *y)  -  0 


(13) 


where 


V  - 


es 


(1+8/2) 


T77 


t  “  8e(l  +  6/2)*^2t 

n 

♦  (1  +  8/2)l/2f±- 


es  1 

y0  "  1,1/2 (Vdi/vi)/(l  +  */2) 

at 

and  the  subscripts  on  $  denote  a  derivative  with  respect  to  that 
variable*  The  quantity  y  represents  the  wave  damping  due  to  electrons, 

*  g 

which  can  result  either  from  resonant  particles  in  a  finite  8  plasma  or 
colllsional  viscosity.  Equation  (13)  is  only  valid  for  Yq  <  1  since  the 
adiabatic  ion  response  can  only  be  Justified  in  this  limit.  Linearizing 
this  equation,  we  obtain  the  complex  eigenvalue  (in  our  normalized  units). 


*  -  k  Ik  |  kZ  Y 

»  ■  «k  ■  -  ^  +  lYn  — *2^7  “  1  — %  * 

^  14k2  0  (14k2)3  1+k2 

Since  Eq.  (13)  is  second  order  in  time,  there  is  an  additional  root. 


(14) 


6 


’X— 


t 


* 


i  *  ilkl(l  +  k2)/Y0,  (15) 

which  Is  spurious  since  it  violates  the  assumption  that  «  kv^  This 
root  Is  growing  and  therefore  must  be  eliminated  before  carrying  out  our 
numerical  computations  (see  Sec.  III). 

In  the  flute  limit  considered  here,  the  magnetic  perturbations  and 

equilibrium  VB  do  not  structurally  alter  the  equation  for  a.  They  enter 
z 

the  equations  through  the  factors  (1  +  g/2)  which  appear  in  the  normalized 
variables  defined  in  Eq.  (13).  We  ignore  this  trivial  finite  0  modifica¬ 
tion  to  the  equations  in  the  remainder  of  this  paper. 

In  the  limit  ♦  0,  Eq.  (13)  reduces  to  the  Hssegawa-Mima 

equation  in  which  the  nonlinearity  arises  from  the  nonlinear  polarization 
drift. ^  This  equation  has  two  invariants,  energy  and  enstrophy,  neither 
of  which  is  preserved  in  the  more  general  Eq.  (13).  When  is  finite,  the 
E  x  B  nonlinearity  also  appears  in  Eq.  (13).  This  nonlinearity  has  been 
considered  previously  in  studies  of  universal  mode  turbulence. 

Our  calculation  differs  from  previous  work  in  that  we  do  not  make  the 
quasi-linear  hypothesis  that 

♦t  +  •  V*  (16> 

in  the  terms  proportional  to  Yq  in  Eq.  (13).  Thus,  our  equation  is  second 
order  in  time  rather  than  first  order.  In  Sec.  Ill  it  is  shown  that  the 

as 

spectrum  always  saturates  when  +  -  1  so  that  the  linear  and  nonlinear  terms 
in  Bq.  (13)  are  comparable.  The  quasi-linear  hypothesis  can  therefore  not 
be  justified  in  drift-wave  turbulence,  k  similar  second  order  equation  has 
been  derived  for  dissipative  drift  waves. 

It  is  recognised  that  both  universal  and  dissipative  drift  wave 
turbulence  are  inherently  three-dimensional  since  k(  «  k*B/E  is  required 
for  the  mode  to  exist  and  the  spatial  variations  in  the  two  directions 
perpendicular  to  B  are  required  to  calculate  the  node  coupling.  In  2-D 
models  of  this  turbulence  the  coupling  of  modes  with  differing  k^s  is 
ignored. 23-26  xhere  is  mo  theoretical  justification  for  this  procedure  so 
that  ears  mast  be  taken  in  applying  the  results  of  model  calculations  to 
experimental  observation «•  By  ce  trust,  the  lower-hybrid-drift  instability 


is  a  flute  mode  so  that  the  neglect  of  the  coupling  along  8  can  be 
justified  and  the  results  of  the  2-D  model  can  be  directly  compared  with 
experimental  observations. 

A  qualitative  understanding  of  the  direction  of  the  flow  of  energy 
described  by  Eq.  (13)  can  be  obtained  by  calculating  the  stability  of  a 


single,  large-amplitude  wave  $q  with  kg  and  u>q  satisfying  the  linear 
dispersion  relation  in  Eq.  (14).  A  perturbation  of  wavevector  k  is  coupled 
through  the  puap  to  nodes  with  k  +  p  k.  (p  -  1,2,...).  For  simplicity,  we 

A  V  "  ^  A 

consider  only  the  coupling  of  (u,  k)  with  its  nearest  neighbors  (u. ,  k, ) 

AAA  A  A  A  V  J  ^ 

where  w+  •  »  +  Ug  and  »  k  +  kg,  i.e.,  p  ■  1.  The  dispersion  relation 

for  this  (k,  k  )  coupled  system  is  (see  Appendix) 


-  -  -  M(iSo»k+)M(k»ko)  - 

*(-.  y  ♦  U0IM  •  °  -«■  fco  *  -k0] 

c(w+,k+) 


where 


e(«,k)  -  i[l  +  kZfl-i^)]  -  ky  +  iye 

\  •  Y0(ky-o)/k3 

M(kl»k2)  “  kl  x  k2*ez [kl  ^  “  i5kl)  ~  k2^  ”  ifik2^ 


When  Yq  and  ye  are  neglected  in  Eq.  (17)  and  we  take  the  limit  a  »  to^, 
the  dispersion  relation  simplifies  to 


“2 

tii 


k2-k2  k2-k2 

u0|2<i> *  i-v2  -V  (-V1 

*  14k  14k  7 


T2  _T2 

+  ko^-  , 

4-  - - - J. 

14k 


(18) 


The  decay  modes  are  purely  growing  with  a  growth  rate  which  peaks  around 

A  A  A  A  A 

k»kg  •  0  with  k/kg  «  0.6  -  0.8,  depending  on  the  magnitude  of  kQ.  A 

necessary  requirement  for  instability  is  that  one  of  the  three  decay  modes 
*  -  27 

(k,  jO  has  a  longer  wavelength  than  the  pump.  This  decay  process  with 
v.  'T  *  . 

k»kg  ■  0  and  k/kg  <  1  is  clearly  seen  in  our  numerical  simulation  during 
the  initial  phase  of  saturation  (Sec.  Ill)  of  the  instability.  For  this 
situation,  in  which  Eq.  (13)  reduces  to  the  Hasegawa-Mima  equation,  the 

wave  energy  inevitably  cascades  to  longer  and  longer  wavelengths  so  that  no 

28 

stationary  wave  spectrum  can  result. 


When  is  included,  this  conclusion  no  longer  remains  valid.  Taking 
the  limit  k^  «  k  and  again  assuming  u>  >>  the  dispersion  relation  is 
given  by 


*2  vs'2;4 ,  v(2+£2i, 

“  — 


Equation  (19)  yields  a  dissipative  instability  when  Yq  is  finite.  The 
destabilizing  term  can  be  traced  back  to  the  E  x  B  nonlinearity.  This 
dissipative  instability  produces  a  flow  of  energy  from  long  to  short 
wavelength  and,  as  will  be  demonstrated,  prevents  the  condensation  of  wave 


energy  at  long  wavelengths,  enabling  the  wave  spectrum  to  reach  a  steady 


state. 


The  dispersion  relation  in  Eq.  (19)  contains  no  amplitude  threshold 
for  instability  as  long  as  k  is  in  the  proper  direction.  In  the  more 
general  dispersion  relation  in  Eq.  (17),  where  the  linear  frequencies  are 
retained,  must  exceed  a  threshold  before  the  instability  onsets  unless 
the  frequency  matching  condition 


\  "  “k+  ±  “0 


is  satisfied.  For  kg  *  k^e  and  k^,  k  «  1,  this  matching  condition  becomes 


“2  *  -  V2  +  Vo)' 


The  locus  of  modes  which  satisfies  Eq.  (20),  those  with  k  <  0  and  Ik  I  < 

•  y  y 


k.,  are  shown  in  Fig.  2.  There  is  no  threshold  for  Instability  of  these 


modes  (neglecting  dissipation).  Nevertheless,  for  the  lower  hybrid  drift 
instability  both  of  the  decay  modes  k  and  k  +  k.  shown  in  Fig.  2  are 
individually  unstable  so  that  the  mode  coupling  process  can  not  saturate 
the  instability. 


Interacting  modes  with  Au  *  0  can  be  driven  unstable  when  exceeds  a 

*****  V 

threshold.  For  the  case  k^  “  k^  e^,  this  threshold  can  be  easily  calcula~ 
ted.  The  dispersion  relation  in  Eq.  (17)  simplifies  to 


-  «V+)2  -  2 


'♦o|2fctek8(kS  -  k2) 


(1  +  k2)(l  +  k2Q+  k2) 


^  V  ‘  ’  V  .  "  V '  .  *  .  ;*•  •*,  .\  »*.  .* 


Fig.  2:  The  locus  of  modes  k  -  kp^g  which  satisfy  the  frequency  matching 

condition  Au»  ■  ■  0,  with  k^  ■  k  +  k^  and 

u.  ■  k  /(I  +  k  ),is  shown. 

*  y 
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V- .0.5 


where  the  dissipative  sources  and  sinks  have  been  neglected.  In  this  limit 

A  A  AAA 

«  ky  ■  0  and  so  that  Am  ■  u>q  -  u^+  is  the  frequency 

mismatch.  To  have  instability, 

A 

Ul2  >  U„l2  S - x - i-  +  k2  x - X - * —  .  (23) 

2(1  +  k2  +k§)(l  +  k§)(kj  -k2) 

A  A 

For  cases  of  interest  k^  -  1  and  UqIc  «  *25.  At  this  amplitude  the 
nonlinearity  overcomes  the  linear  frequency  mismatch  of  the  modes  and  a 
broad  range  of  modes  is  driven  unstable.  As  a  consequence,  we  would  expect 
the  wave  spectrum  to  saturate  at  around  this  level  even  when  «  1,  i.e., 
the  driving  rate  is  small.  This  conclusion  is  verified  in  our  numerical 
computations . 

The  particle  flux  can  be  calculated  self-consistently  with  the  wave 
spectrum.  The  rate  of  change  of  the  density  n(x)  =  <ne(x,y)>y,  where  <  >y 
is  a  spatial  average  over  y,  is  obtained  by  averaging  the  electron 
continuity  equation, 

|S.  +  |_<nV  >  -0.  (24) 

3t  3x  e  ex  y  v 

To  lowest  order,  *  -  (c/B)3$/3y,  the  polarization  drift  being  small. 
The  electron  density  can  not  be  calculated  directly  since  Eq.  (24)  for  ne 
can  not  be  inverted  analytically.  However,  since  ne  =  n^,  we  can  use  the 
expression  for  n^  in  Eq.  (4)  to  obtain  an  explicit  expression  for  the 
diffusion  coefficient  D  -  -  <neVex>y/  (3n/3x), 

£  ■  m0  -  (f)1/2  *  4_  - 1  (25) 

D0  ■  °.s“lh 

Note  that  the  normalization  factor  Dq  of  the  diffusion  coefficient 
corresponds  to  a  step  size  p  with  a  correlation  time  oj,\  In  the 

a  a  in 

quasilinear  limit  3/aT  ■  -  iu^,  and  D  is  given  by 

V  "  (f)1/2  I  k  kJUJ2*1  +  ^2>_1*  <26> 

A  A 

The  time  evolution  of  both  D  and  will  be  computed  along  with  the  wave 
spectrum  in  Sec.  III. 
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III.  Numerical  Computations 


The  nonlinear  equation  for  4  in  Eq.  (13)  is  too  complex  to  solve 
analytically  and  we  have  therefore  adopted  a  computational  approach.  It  is 
convenient  to  separate  the  second  order  Eq.  (13)  into  two  complex  first 
order  equations  by  defining  a  new  function 

♦  iW  +  Y0Iv“1I(Jt  +  %y).  (27) 

A 

Equation  (13)  for  4  can  then  be  written  as 

A  A  AA  A  A 

♦  "*+♦  +  Y.4>  +  V*  x  e  .V*  -  0.  (28) 

I  T  7  “  Z 

Since  Yq  «  1,  we  find  from  Eq.  (27)  that  to  lowest  order 

♦  *  V~2*.  (29) 

The  second  term  on  the  right  side  of  Eq.  (27)  is  already  small  since  it  is 

A 

proportional  to  y0  80  that  $  in  this  term  can  be  approximated  to  lowest 
order.  Equation  (27)  then  becomes  an  evolution  equation  for 

+T  -  -  +y  +  Ivlv2(*  -  V2j)/Y0*  (30) 

Finally*  eliminating  9  in  Eq.  (28)  using  Eq.  (30),  we  obtain  a  similar 

T  A 

evolution  equation  for  $* 

♦T  -  -  iy  +  Ye*  -  1>y  +  Ivlv2(*  -  V2J)/y0  -  V J  x  e8.V+.  (31) 

Note  that  the  iteration  procedure  which  was  carried  out  on  Eq.  (27) 
eliminated  the  spurious  unstable  root.  The  dispersion  relation  which 
results  from  linearising  Eqs.  (30)  and  (31)  yields  the  correct  eigenvalue 
of  the  lower  hybrid  drift  instability  given  in  Eq.  (14)  plus  a  spurious 
damped  root, 


«  -  -  ilk|3<l  +  k2)/Y(),  (32) 

which  docs  not  cause  any  numerical  difficulties.  All  of  our  numerical 

A 

computations  ere  based  on  the  two  coupled  equations  for  <j>  and  in  (30)  and 
(31). 

A 

The  inclusion  of  the  electron  damping  y^  in  Eq.  (31)  is  essential  in 
order  to  obtain  a  saturated  spectrum  since  in  the  absence  of  this  term  all 
modes  in  k  space  are  unstable  (see  the  growth  rate  in  Eq.  (14)].  We  have 
used  two  forms  of  electron  dissipation  in  our  calculations.  The  first 
arises  from  the  VB  resonant  electrons  in  a  finite  g  plasma.  Explicit 
expressions  for  this  damping  rate  have  been  derived  in  Ref.  9.  In  the  long 
wavelength  limit 


%  ■  v>lky|k' 


with 


U  +  28i)2 
8^(2  +  6l) 


exp(-  2/8  ) 


where  6  and  g.  are  the  Individual  electron  and  ion  values  of  beta. 

®  *  A  A  A 

short  wavelength  limit,  scales  as  k^k2  so  we  take  as  a  model 


In  the 


\  “  Wvk4/(l  +  k2)*  (35) 

A 

An  important  feature  of  the  SB  resonant  damping  is  that  it  scales  as  k5  as 

A  A 

k  <*•  0.  Since  the  growth  term  in  Eq.  (14)  scales  as  k3  in  the  long  wave** 

A  A 

length  limit,  modes  with  k  ■  0  remain  unstable  even  as  k  ♦  0.  For  large 
*  ^  . 
k  all  modes  are  stable.  The  boundary  between  stable  and  unstable  regions 

of  k  space  is  shown  in  Fig.  3a  for  Y^q/Yq  "  an<*  1*°*  The  area  below 

the  bounding  curves  is  unstable.  As  Y^q/Yq  increases  the  unstable  spectrum 

collapses  toward  long  wavelength  but  unstable  modes  always  remain. 

Electron  collislonal  viscosity  is  important  in  laboratory  experiments 

where  measurements  of  lower-hybrid-drift  wave  turbulence  have  been  made. 

21 

The  electron  damping  rate  for  this  case  is  given  by** 


irk  '-m-  ^  ^  fin'  ^  ■■{Jr<  ^  <.i<^  -.-  * 


Y 


e 


where 


veVnep 


es’ 


(36) 


(37) 


v  is  Che  electron-don  collision  frequency 
© 

v  -  2.9  x  l<f6(Xn  /T  3/2)s'\ 

“  C  © 

X  is  the  Coulowb  logarithm,  n#  is  the  electron  density  in  ca  and  Te  is 
the  electron  teaperature  in  ev.  In  this  case  the  collisional  damping 

A 

dominates  the  growth  at  both  large  and  saall  k  ao  that  the  unstable  region 

as  as 

of  k  space  is  localised  around  k  «  1  aa  shown  in  Fig.  3b  for  several  values 
of  Vq/yq.  There  is  no  Instability  when  vq/yq  >  /3/16. 

as 

The  coupled  equations  for  4  and  4  in  (30)  and  (31)  are  solved  using  a 
pseudo-spectral  method  code  developed  by  Fyfe  et  al. ,  based  on  the  work 
of  Or a sag. 30  The  dependent  variables  4  and  4  are  Fourier-decomposed , 

4(x,t)  -  I  4(k , t )exp( ik. x) ,  (38) 

k 

ma 

A  as  a  * 

where  k  ■  n/X  and  n  »  n  e  +  n  a  with  1^  and  n_  integers.  The  parameter  X 
fixes  the  value  of  |n|  for  which  Ik]  •  1.  The  nonlinear  term  in  Eq.  (31) 

V  V  AA  A  A 

is  computed  by  fast  Fourier  transforming  74  and  V4  from  k  to  configuration 

asA  A  A  V 

apace,  calculating  the  product  74  x  a  *74  in  configuration  space  and  then 

•  A 

faat  Fourier  transforming  the  result  back  to  k  space.  The  equations  are 
then  stepped  forward  in  time  in  k  space  and  the  cycle  is  repeated. 

The  numerical  results  presented  in  this  paper  are  nominally  computed 
on  a  32  x  32  mesh.  However,  to  prevent  aliasing  of  the  wave  energy  during 
the  fast  Fourier  transformation,  it  is  nscessary  to  sero  all  modes  with  n  > 
32/3  so  that  the  useful  voltse  of  k  space  is  actually  much  less  than  one 
would  expect. 

A  number  of  tests  have  bean  made  to  ensure  that  the  code  is  correctly 
advancing  the  equations  in  time.  The  growth  rates  of  the  nodes  as  obtained 
from  the  code  during  the  linear  growth  phase  have  been  checked  with  the 


,T.  VWV'  « 
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Th«  region  of  instability  in  k  space  (below  bounding  curves)  is 
shown  for  (a)  collisionless  plasaa  and  (b)  colllsional  plains, 
where  Yq»  Y#q  *nd  vo  *re  th®  nomalised  ion  growth  rate, 
electron  VB  danplng  rate  and  electron  colllsional  damping  rate. 
Colllsional  dasping  stabilizes  the  spectrua  for  Y^/Yq  >  /T/16. 


solution  of  the  linear  dispersion  relation.  We  have  zeroed  all  but  three 
interacting  nodes  in  the  full  code  and  cross-checked  their  nonlinear 
behavior  with  an  independent  code  written  to  specifically  study  the 
interaction  of  only  three  waves.  The  sensitivity  of  our  results  to  the 
step  size  of  the  tine  integration,  as  well  as  to  the  number  of  modes  (some 
runs  have  been  carried  out  on  a  64  x  64  mesh),  have  also  been  checked. 

IV.  Numerical  Results:  Collisionless  Damping 

In  collisionless,  finite-g  plasmas  the  VB  resonant  electron  damping 
limits  the  range  of  unstable  modes  as  shown  in  Fig.  3a.  In  this  section  we 
present  the  results  of  our  numerical  computations  for  this  case.  Pre¬ 
liminary  results  were  previously  presented  for  this  collisionless  case 
although  the  structure  of  our  electron  damping  is  now  more  realistic  than 

Oft 

that  which  was  used  in  past  computations . 

In  Figs.  4-6  we  present  in  some  detail  the  results  of  our  computations 

A 

for  the  case  -yn  ■  0.5,  -y  n  ■  0.01  and  X  ■  5  (k  p  “1  for  n  •  5).  The  4 
spectrum  is  initialized  with  random  noise  with  10  3  and  Eqs.  (30)  and 
(31)  are  evolved  until  the  wave  energy  (in  our  normalized  units), 

5  -i  £  (1  +  k*)|l|2,  (39) 

k 

A 

approaches  a  steady  state  value.  The  time  history  of  W,  the  root-mean- 
square  potential, 

P  -  <}2>1/2  =  (Ii;k|2)1/2  (40) 

k 

A  A 

and  the  diffusion  coefficients  D  and  [Eqs.  (25)  and  (26),  respectively] 
are  shown  in  Fig.  4.  All  of  these  quantities  exhibit  a  similar  temporal 
behavior.  The  initial  decay  (r  <  4)  is  associated  with  the  rapid 
dissipation  of  energy  initialized  in  the  damped  modes,  and  is  followed  by  a 
linear  growth  phase  (4  <  t  <  30).  Subsequently,  mode  coupling  occurs  which 
leads  to  saturation  of  the  instability,  albeit  with  some  initial 
overshoot  (t  -  30).  The  levels  of  the  total  wave  energy  and  other 
parameters  of  Fig.  4  are  quite  stationary  in  time  after  saturation.  Also, 
the  stationary  values  of  all  four  quantities  are  relatively  Insensitive  to 

A 

the  initialization  of  4. 


luivunnrv.  r.  ^LT7.'7 ■*  ■?  vj  r-  «r  '  ■■  7 


'.•  *■■* t'.*  o  •:*  -'.*  v* •: 


In  Fig.  4  the  potential  P  asymptotes  to  2.0  after  saturation  so  that 

* 

4  -  1,  which  is  consistent  with  the  discussion  in  Sec.  II.  The  quasilinear 

*  A 

diffusion  coefficient  (D  , )  tracks  the  actual  diffusion  coefficient  (D) 
quite  well  during  the  entire  time  evolution  of  the  instability.  An 

A 

important  point  which  must  be  emphasized  with  regard  to  D  is  that  both 
species,  electrons  and  ions,  continue  to  exchange  both  energy  and  momentum 
even  after  a  steady  state  is  reached;  the  electrons  through  the  VB  reso¬ 
nance  and  the  ions  by  direct  resonant  interaction.  If  the  instability  had 
saturated  by  ion  trapping  and  the  electrons  had  no  resonant  interaction 
with  the  wave,  there  would  be  no  diffusion  in  the  steady  state  since  the 
electrons  could  not  exchange  momentum  with  the  ions.  Both  species  must 
have  a  dissipative  interaction  with  the  waves  to  have  diffusion. 

In  Fig.  5  we  show  a  sequence  of  snapshots  of  the  2-D  wave  spectrum  in 
n  space  as  the  instability  grows  and  saturates.  The  times  at  which  the 
snapshots  are  taven  are  indicated  by  the  arrows  in  Fig.  4a.  Only  ^  >  0  is 
shown  since  the  spectrum  for  n^  <  0  can  be  obtained  from  nx  >  0  by  the 
reality  condition  $(-  n)  ■  +(n).  Figure  5a  shows  the  spectrum  during  the 
linear  phase  of  the  instability.  The  spectrum  is  strongly  peaked  around 
the  most  unstable  modes,  (n^,  ny)  -  (0,  +  5).  Figure  5b  shows  the  spectrum 
just  prior  to  saturation.  Note  the  development  of  the  prominent  peak  at 
(4,  0),  which  is  a  marginally  stable  mode  since  riy  ■  0,  and  the  broadening 
of  the  main  peaks  of  the  spectrum  in  the  ^  direction.  These  results  are 
entirely  consistent  with  the  discussion  of  the  parametric  excitation  of 
daughter  waves  by  a  pump  in  Sec.  II.  The  pump  wave  (0,  +  5)  couples  and 
destabilizes  the  (4,  0)  and  (4,  +  5)  modes.  These  secondary  peaks  quickly 
grow  to  large  amplitude  and  excite  other  modes  in  turn.  This  mode  coupling 
process  culminates  as  the  total  wave  energy  saturates  and  spreads  through 
most  of  the  unstable  or  weakly  damped  volume  of  n  space  as  shown  in  Fig. 
5c.  Two  significant  features  of  Fig.  5c  are:  the  shift  in  the  peaks  of 
the  wave  spectrum  toward  long  wavelength  (0,  +2);  and  the  nearly  isotropic 
spectrum  of  waves  surrounding  these  peaks.  The  detailed  wave  spectrum 
exhibits  substantial  variability  in  time,  even  after  saturation  when  the 
total  wave  energy  is  nearly  constant.  This  can  be  clearly  seen  by 
comparing  the  spectrum  in  Fig.  5c  with  that  at  a  later  time  in  Fig.  5d. 
The  unstable  and  stable  modes  clearly  continue  to  exchange  energy  in  a 
dynamic  manner  even  after  saturation.  Nevertheless,  the  wave  spectra  at 
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(c)  (d) 

5:  Snapshots  of  the  two-dimensional  wave  spectra  are  shown  (a) 

during  the  linear  phase,  (b)  during  saturation,  (c)  just  after 

saturation  and  (d)  well  after  saturation  for  the  run  shown  in 

* 

Fig.  4,  where  a  »  5  k.  The  times  of  the  snapshots  are  marked 
with  arrows  in  Fig.  4a. 
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late  time  are  always  characterised  by  peaks  which  are  shifted  towards  long 
wavelength  from  the  most  unstable  mode  and  by  a  broad  spectrum  of  noise . 
Finally,  in  Fig.  6a  we  show  the  wave  spectrum  averaged  over  the  time 
interval  112  <  t  <  192.  Much  of  the  irregularity  which  was  evident  in  the 
instantaneous  spectra  has  now  disappeared.  The  averaged  spectrum  is 
strongly  peaked  around  (0,  +  3)  with  much  smaller  secondary  peaks 
at  (+  2,  0).  A  contour  plot  of  this  average  spectrum  is  shown  in  Fig. 
6b.  Note  the  cleft  in  the  spectrum  along  ky  -  0. 

A  number  of  runs  have  been  made  with  different  values  of  the  drift 
parameter  Yq  and  the  coefficient  of  the  electron  damping  Yeo*  *n  F*8*  ?  we 
show  the  diffusion  coefficient  D  and  potential  P  at  saturation  as  functions 
of  the  drift  parameter  Yq*  The  electron  damping  y6q  was  varied  with  y0  80 
that  the  spectrum  of  unstable  waves  in  k  space  did  not  change,  l.e.,  y6q  “ 
0.1  Yq*  The  potential  P  is  quite  insensitive  to  Yq  while  the  diffusion 
coefficient  scales  approximately  as  Yq  (tl»*  reference  curve  D  «  0.76  Yq  is 
shown  for  comparison).  The  insensitivity  of  P,  and  consequently 

A 

4> ,  to  Yq  1*  consistent  with  the  idea,  expressed  previously  in  Sec.  II,  that 
the  non-linear  polarization  drift  in  Eq.  (13)  must  exceed  the  linear 
frequency  mismatch  of  the  Interacting  waves  in  order  for  effective  energy 
exchange  (and  therefore  saturation  of  the  instability)  to  take  place. 

A  A 

In  Fig.  8  we  show  P  and  D  as  a  function  of  yq  with  y6q  ■  0.11.  In 
this  case  the  linearly  unstable  spectrum  collapses  towards  longer  wave¬ 
lengths  as  Yn  decreases  as  previously  shown  in  Fig.  3a.  The  dependence  of 

A  A  " 

P  and  D  on  Yq  in  this  case  is  basically  similar  to  that  in  Fig.  7.  The 
normalized  potential  P  has  a  weak  dependence  on  Yq*  increasing  slightly  as 
Yq  decreases.  The  diffusion  coefficient  again  scales  as  y§  (the  reference 
curve  D  ■  0.76  y§  is  shown). 

The  somewhat  surprising  conclusion  which  can  be  drawn  from  Figs.  7  and 

A 

8  is  that  the  diffusion  coefficient  D  is  very  Insensitive  to  the  dissi- 

A 

pat  ion  rate  y€q5  much  less  sensitive,  for  example,  than  P.  For  Yq  “  0.25, 
the  diffusion  coefficient  D  -0.06  in  both  Figs.  7  and  8  even  though  ycq  ■ 
0.03  and  0.11,  respectively.  By  contrast,  P  is  50%  higher  in  Fig.  8  than 

A 

in  Fig.  7  for  Yq  “  *25.  As  ycq  1*  increased  for  a  fixed  value  of  Yq»  4> 
increases  while  the  spectrum  shifts  towards  longer  wavelength  so  that  D 

A 

does  not  change  [see  Eq.  (26)  for  Dql].  Th«  shift  towards  long  wavelength 


Fig.  7;  The  dependence  of  the  diffusion  coefficient  D  and  the  root-mean- 

A 

square  potential  P  on  the  ion  driving  term  Yq  *  /w  v<ji^vi  is 
shown  for  a  fixed  linearly  unstable  spectrum  of  waves 
(YeQ/T0  “  0.1  and  Vq  ■  0).  The  reference  curve  is  given 
by  D  -  0.76  y02- 


A  A 

Fig.  8:  The  dependence  of  D  and  P  on  Yq  is  shown  for  a  fixed  electron  VB 


daaping  rate 


TJ0 


■  0.11  with  v. 


0.  The  reference  curve  is 


again  given  by  D  ■  0.76  yq5 


28 


as  YeQ  increases  can  be  clearly  seen  in  plots  of  the  2-D  wave  spectrum.  By 
comparing  the  quasilinear  diffusion  coefficient  in  Eq.  (26)  with  the 
portion  of  the  linear  growth  rate  in  Eq.  (14)  which  comes  from  the  resonant 
ion  interaction  and  the  expression  for  the  wave  energy  in  Eq.  (39),  we  find 
that  roughly 


D  -  8Wi/3t, 


where  3W^/dt  is  the  rate  at  which  energy  is  extracted  from  the  ions.  Thus, 
our  numerical  results  indicate  that  the  rate  at  which  the  instability  is 
taking  energy  from  the  ions  is  nearly  independent  of  the  electron 
damping.  Of  course,  damped  modes  must  always  be  included  in  the  simulation 
or  the  wave  energy  does  not  saturate.  The  damping  rate  y6q  “ust  therefore 
remain  finite. 

V.  Numerical  Results:  Collisional  Damping 

In  laboratory  and  ionospheric  plasmas  classical  damping  can  dominate 
the  collisionless  VB  electron  damping.  As  noted  previously  in  Fig.  3b, 
collisional  damping  stabilizes  both  long  and  short  wavelength  modes  and, 
unlike  the  VB  resonant  electron  damping,  can  stabilize  the  ent ire 
spectrum.  The  results  of  our  numerical  computations  for  this  ca*£  a.e 
summarized  in  Fig.  9  where  we  show  the  dependence  of  the  saturated  values 

A  A 

of  D  and  P  on  the  collisional  damping  coefficient  Vq  for  Yq  ■  1.0.  The 
diffusion  coefficient  D  and  potential  P  are  virtually  independent  of  Vq. 
This  result  is  rather  surprising  because  as  Vq  increases  the  range  of 
unstable  modes  in  k  space  decreases  dramatically  (Fig.  3b)  and,  as  a 
consequence,  the  characteristic  wave  spectra  after  saturation  in  the  three 
cases  presented  in  Fig.  9  are  rather  different.  The  spectra  in  the 
cases  Vq  ■  0.05  and  0.10  are  similar  to  those  presented  in  Sec.  IV,  the 
spectrum  being  somewhat  broader  for  Vq  ■  0.05  because  the  range  of  linearly 

unstable  modes  extends  to  larger  values  of  Ik  I .  For  vQ  ■  0.15  the  linearly 

unstable  region  of  k  space  has  become  rather  localized  and,  as  shown  in 

Fig.  10,  the  saturated  wave  spectrum  is  much  more  peaked.  The  two  dominant 

peaks  at  n  -  (0,  +  4)  correspond  to  the  most  unstable  modes  while  the 
secondary  peak  at  n  •  (2,0)  is  the  parametrically  generated  daughter 
wave.  In  this  case  the  damping  of  the  modes  with  finite  is  so  large 


"Hint;; 


(except  for  tiy  •  0)  that  the  energy  which  Is  pumped  into  them  is  rapidly 
dissipated  and  the  wave  spectrum  remains  quite  narrow.  Nevertheless,  the 

A  A 

saturated  values  of  both  P  and  D  are  virtually  the  same  as  those  computed 
for  smaller  values  of  Vq«  This  result  is  a  consequence  of  the  structure  of 
the  nonlinearity  in  Eq.  (13)  which  only  couples  modes  k",  k  with 
k"  x  k  *  0.  For  the  case  of  a  narrow  wave  spectrum  such  as  that  shown  in 

v  v 

Fig.  10  the  coupling  of  the  modes  is  quite  weak.  Therefore,  even  though 
the  growth  rates  of  the  unstable  waves  are  small,  the  coupling  of  wave 
energy  to  the  stable  modes  is  correspondingly  very  weak  and  the  saturation 
amplitude  remains  large. 

From  our  numerical  results  we  '  conclude  that  the  saturation  of  the 
lower-hybrid-drift  instability  remains  large  even  as  collisions  force  the 

A 

mode  towards  marginal  stability.  At  marginal  stability,  of  course,  4  0. 

This  rather  discontinuous  behavior  as  collislonality  increases  seems  rather 
unphyslcal  and  compels  us  to  re-evaluate  our  assumptions.  In  the 

Introduction  we  argued  that  ion  trapping  was  not  a  viable  saturation 
mechanism  when  a  broad,  two-dimensional  wave  spectrum  is  excited.  As  we 
approach  marginal  stability  by  increasing  the  collislonality,  the  spectrum 
becomes  quite  narrow  so  that  ion  trapping  should,  at  some  point,  become 
effective  and  saturate  the  instability  at  a  lower  amplitude  than  that  which 
we  have  calculated  from  our  mode  coupling  theory.  Ion  trapping  occurs  when 
the  bounce  time  u.  .  of  an  ion  in  the  wave  potential  becomes  comparable  to 
the  growth  time  y”^of  the  Instability  or 


wJi  “  •  Y2 


In  our  normalized  units. 


♦  ■  Y2(«1/%)1/2Vdl/v1 


for  k  -  1  so  that  $  ♦  0  as  the  mode  approaches  marginal  stability. 
Saturation  by  ion  trapping  occurs  at  e  lower  amplitude  than  that  resulting 
from  mode  coupling  when 


T  <  C,/»1)1/4(V1/V<)1)1/2. 
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VI.  Sunary  and  Applications 


A  nonlinear  wave  equation  [see  Eq.  (13)]  has  been  derived  which 
describes  two-dimensional  (k»B  •  0)  lower-hybrid-drift  wave  turbulence  in  a 
finite  6  plasaa.  As  in  the  linear  theory  of  this  instability  in  the  weak 
drift  regime,  the  ions  are  treated  as  unmagnetized  and  to  lowest  order 
respond  adiabatically  to  the  potential  fluctuations.  The  resonant  ion 
contribution  is  included  as  a  correction  to  provide  the  driving  energy  for 
the  instability.  The  nonlinear  portion  of  the  wave  equation  arises  from 
the  electrons  through  their  E  x  and  polarization  drifts.  In  addition, 
electron  dissipation  is  included  either  through  the  collisionless  VB  reso¬ 
nance  or  through  colllsional  viscosity. 

The  nonlinear  wave  equation  has  been  solved  numerically  using  a 
pseudo-spectral  method  code  to  obtain  the  evolution  and  saturation  of  the 
wave  spectrum  as  well  as  the  self-consistent  particle  flux.  The  wave 
energy  typically  saturates  as  the  peaked  unstable  spectrum  spreads 
throughout  the  unstable  and  weakly  damped  regions  of  k  space.  The 
saturated  wave  spectrum  is  characterized  by  two  main  components:  a  broad, 
nearly  Isotropic  spectrum  with  |k|  >  p  !;  and  a  jagged  spectrum  consisting 
of  narrow  peaks  in  the  range  |k|  -  0.5  p  .  These  peaks  are  typically 
displaced  towards  long  wavelength  from  the  linearly  most  unstable  waves  by 
almost  a  factor  of  two.  The  peaks  move  around  in  k  space  as  the  waves 
exchange  energy  in  a  dynamic  fashion. 

Our  numerical  computations  demonstrate  that  the  wave  energy  saturates 
when  the  nonlinearity  becomes  comparable  to  the  linear  frequency  of  the 

modes , 


dt 


01  ■«*  V  *V  —  kV  , 
Tc  ^e  e  ' 


(44) 


where  is  the  E  x  B  velocity  of  the  electrons  and  ■  lcy ( 1  + 
k2p2#).  For  the  most  unstable  modes,  kp#>  -  1,  -  Vai/Peg  ~  “ihVdi^vi 

and  (44)  can  be  written  as 


(45) 


et/T±  »  2-3  (2ae/oi)1/2Vdl/v1  , 

where  the  numerical  factor  comes  from  the  detailed  computations.  We  have 
also  found  that  the  saturation  level  is  relatively  Insensitive  to  both  the 
magnitude  of  the  electron  dissipation,  as  long  as  the  damping  is  sufficient 
to  stabilize  the  shortest  wavelengths  in  the  computational  grid,  and  the 
damping  (growth)  rate  of  the  long  wavelength  modes. 

Our  interpretation  of  the  saturation  amplitude  in  Eq.  (45)  is  that 
effective  energy  exchange  between  the  linearly  growing  and  damped  modes  can 
only  take  place  when  the  nonlinearity  becomes  sufficiently  large  to 
overcome  the  linear  frequency  mismatch 

"  “kl  +  “k2  +  “k3  (46) 

between  Interacting  waves,  k^,  and  kj.  Below  the  threshold  given  in  Eq. 
(45),  only  a  few  modes  which  satisfy  Au  *  0  can  exchange  energy.  Above 
this  threshold  all  modes  can  exchange  energy.  The  saturation  amplitude  in 
Eq.  (45)  also  corresponds  to  the  threshold  for  electron  E  x  B  trapping-’ ^  so 
that  strong  mode  coupling  is  really  a  consequence  of  the  strong  nonlinear 
behavior  of  the  electrons  during  this  trapping  process. 

The  particle  flux  has  been  confuted  self-consistently  in  parallel  with 
the  nonlinear  evolution  of  the  wave  spectrum.  Both  the  exact  and 
quasillnear  expressions  for  the  particle  diffusion  coefficient  D  are  then 
calculated  from  the  flux.  The  quasillnear  expression  (D^)  tracks  the 
exact  expression  (D)  within  a  factor  of  two.  To  lowest  order 
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where  the  numerical  factor  cones  from  the  detailed  computations.  The 
diffusion  coefficient  in  (48)  is  very  insensitive  to  the  electron 
dissipation,  less  sensitive  than  the  saturation  amplitude  e^/T^  given  in 
Eq.  (45).  At  first  glance  Eq.  (47)  seems  to  imply  that  the  usual  estimate 
Y  ■»  k2D  or  D  .  p2  y,  where  y  is  the  maximum  growth  rate  of  the  instability, 
is  quite  accurate.  However,  the  insensitivity  of  D  to  the  electron 
dissipation  implies  that  this  estimate  is  not  valid.  The  growth  rate  y 
decreases  rather  sharply  as  the  electron  dissipation  is  increased  so  that 

D  -  o2  y  should  also  decrease.  It  does  not  so  that  this  estimate  falls  to 

es 

reproduce  the  computational  results. 

We  now  compare  the  results  of  our  mode  coupling  calculation  with 
laboratory  observations  of  lower-hybrid-drif t  wave  turbulence  and  with  the 
results  of  recent  computer  simulations  of  this  instability.  A  CO2  laser 
scattering  experiment  was  recently  carried  out  on  the  Garching,  10m 
0-pinch.^  The  measured  density  fluctuations  were  flute-like  (k«B  *  0) 
with  clear  peaks  around  k_  -  0.5  p  1.  The  lower-hybrid-drift  instability 

J  ®8 

was  therefore  identified  as  the  source  of  these  fluctuations.  Our  computa¬ 
tional  results  are  in  good  agreement  with  several  prominent  features  of  the 
observed  wave  spectra.  They  measure  a  broad  angular  spectrum  with  k^  <  lty. 
They  also  observe  that  the  peak  in  the  wave  spectrum  is  shifted  towards 
long  wavelength  by  about  a  factor  of  two  from  the  linearly  most  unstable 
mode.  In  the  experiment  data  were  taken  for  three  different  filling 
pressures,  corresponding  to  three  values  of  W*i-  We  have  found  that 
classical  electron-ion  collisions  have  a  significant  influence  on  the 
growth  rate  of  the  lower-hybrid-drift  Instability  in  this  experiment 
especially  for  the  two  highest  filling  pressures.  In  Table  I  we  summarize 
some  of  the  parameters  of  their  experiment  Including  the  ratio  of  the 
electron-ion  collision  frequency  to  the  critical  collision  frequency, 

v.  •  -8U(V’i)2  “lh 

required  to  completely  stabilize  the  lower-hybrld^drift  wave.  In  calculat¬ 
ing  vai  we  have  Ignored  impurities  by  taking  Zeff  -  1.  The  highest 
pressure  case  (8m  Torr)  is  very  close  to  marginal  stability  and,  given  the 
uncertainty  in  experimental  parameters,  could  even  be  stable.  We  have 
completed  numerical  computations  for  the  three  cases  given  in  Table  I.  For 


l  v  iv  rv^  w  .!.■»  ^  '.v  w  .a  a-  ia  \a  i> 


Table  I 

Parameters  of  Garching  Experiment 


P0(m  Torr) 

ne(1014/cm3) 

Te(ev) 

B(T) 

Wvi 

ve(106/sec) 

v  /v 
e  c 

3 

1 

115 

.79 

.52 

3.3 

.0067 

5 

2.3 

74 

.77 

.25 

13.7 

.12 

8 

5.2 

77 

.73 

.17 

28.5 

.56 

the  tvo  lowest  filling  pressures,  3  and  5  m  Torr,  the  instability  saturates 

with  n/n  “  3.4  x  10*~^  and  1.4  x  10"2,  respectively,  compared  with  the 

measured  values  of  1.4  x  10  and  2.3  x  10  .  Our  calculated  saturation 

levels  are  2-3  times  larger  than  measured  experimentally.  In  the  case  of  8 

m  Torr  the  spectrum  is  close  to  marginal  stability  and  as  a  consequence  the 

saturated  wave  spectrum  is  strongly  peaked  around  the  narrow  band  of 

unstable  modes  at  k  p  ~  0.58.  Because  the  spectrum  is  quite  narrow,  ion 

y  es 

trapping  cannot  be  neglected  so  that  the  supposition  made  in  Ref.  5  that 
ion  trapping  saturates  the  instability  is  probably  valid  for  this  case. 

A  number  of  computer  simulations  of  the  lower-hybrid-drift  Instability 
have  been  carried  out  over  the  past  several  years.1®*13,32”3*  Before 
comparing  our  results  with  these  simulations,  we  would  like  to  point  out 
several  difficulties  which  must  be  overcome  in  simulating  this  instability 
which  are  not  widely  recognised.  The  large  disparity  between  the  electron 

gyro  time  scale,  n  ,  and  the  growth  time  of  the  instability,  (v./Vj,.)* 

1/2  —1  *  1  01 
(mi/me)i/,iQe  ,  have  forced  the  theorists  to  use  artificially  small  mass 

ratios,  (s^/m^1/2,  and  high  drift  velocities,  V^/v^  in  their  codes.  As 

a  consequence,  the  parameter 

pJl. ' 

which  should  be  a  small  number,  is  of  order  unity.  Cases  with  larger  mass 

ratios  had  to  be  run  with  larger  drift  velocities  so  Peg/Ln  has  not  been 

varied  significantly.  The  parameter  p  /L  is  significant  for  our 

computations  for  two  reasons.  First,  the  ratio  of  the  rate  of  change  of 

the  magnetic  free  energy  to  particle  drift  energy  in  a  finite  £  plasma 

scales  as  l£/p2  ,13  For  realistic  parameters  the  magnetic  free  energy 

therefore  greatly  exceeds  the  drift  free  energy  and  the  lower-hybrid-drift 

instability  has  no  free  energy  bound.  In  the  simulations  this  has  not  been 

the  case.  Second,  in  calculations  of  the  nonlocal  structure  of  the  lower 

hybrid  drift  instability,  it  has  been  shown  that  the  number  of  unstable 

harmonics  in  the  x  direction  scales  as  L  /p  so  that  for  realistic 

n  **  is 

parameters  a  broad  spectrum  of  unstable  modes  should  be  excited.'9'7  This  is 
again  not  the  case  in  many  simulations.  Recently,  fluid-particle  hybrid 

codas  have  been  developed  and  simulations  of  the  lower  hybrid  drift 

13 

instability  have  been  carried  out  with  p  /L  «  1.  In  these  simulations 

es  n 
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a  broad  spectrum  of  nodes  in  k  space  is  excited  with  k^  <  ky. 
the  wave  spectrum  along  ky  -  0,  as  shown  in  Fig.  6,  is  also  seen  in  the 
wave  spectra  from  these  simulations.  Finally,  the  amplitude  of  Che 
potential  fluctuations  at  saturation  as  obtained  from  these  simulations  is 
given  by 

.,/Tj  .  2.4(2me/.1)1/2(Vdl/»1)J 

which  is  in  excellent  agreement  with  our  results  [see  Eq.  (45)].  One 
caveat  to  this  comparison  is  that  most  of  their  simulations  were  carried 
out  for  Vjj/v^  >  1  while  our  theory  strictly  applies  only  for  W’i  <  >■ 

Finally,  we  now  consider  two  physical  systems,  the  earth's 
magnetosphere  and  ionosphere,  where  lower-hybrid-drift  turbulence  is 
expected  to  be  important.  The  plasma  in  the  neutral  sheet  of  magnetotail 
is  essentially  collisionless .  The  lower-hybrid-drift  Instability  has 

previously  been  proposed  as  a  mechanism  for  dissipating  magnetic  energy  in 

/ 

this  reversed  field  configuration.  In  the  most  simple  l-D  model  of 
magnetic  field  annihilation  the  flux  merging  velocity  V£  is  simply  given 


where  v  is  the  collision  frequency  (anomalous  or  classical)  and  X  is  the 
scale  length  of  the  magnetic  field.  The  anomalous  collision  frequency  due 
to  lower-hybrid-drift  turbulence  is  given  by  v  •  D/p2  so  that  from  Eq. 
(48),  we  find 

V\  -  2-*(«./-1)l/2(V’,i)3> 

where  V2  •  B2  Airman  is  the  Alfven  velocity.  The  rate  of  dissipation  of 
magnetic  flux  is  quite  small  during  typical  quiet  conditions  when  w 
vt  -0.1  while  during  substorm  activity  when  substantial  thinning  of  the 
sheet  has  been  reported  (V^/v^  <  1),  the  dissipation  rate  can  be  quite 

large. 


The  lower-hybrid-drift  instability  has  also  been  suggested  as  a 

Mechanism  to  generate  small-scale  irregularities  (X  <  1  m)  in  the 

ionosphere  during  equatorial  spread  F  (ESF),^»*^  and  observational 
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evidence  supports  this  mechanism*  ’  An  important  quantity  to  expert- 

O  Q  IQ 

mentalists  is  the  power  spectrum  of  density  fluctuations  during  ESF. 
Recently,  Singh  and  Szuszczewicz  have  presented  a  composite  of  ESF  wave 
spectra  from  medium  to  short  wavelengths.  They  find  that  (1)  the  spectra 
in  the  medium  and  intermediate  wavelength  domain  (X  ~  50  km  -*•  200  m)  scale 
as  k  and  k  respectively;  (2)  the  transitional  wave¬ 
lengths  (X  ~  200  m  ♦  20  m)  have  a  k  dependence,  presumably  due  to 

drift  waves;  (3)  the  power  spectrum  breaks  at  kp^  >  1  (x  <  20  m)  to  a 
shallower  k  dependence;  and  (4)  the  shortest  wavelength  observations 
(1  m  <  x  <  20  m)  are  consistent  with  a  wave-particle  interaction  such  as 

v  v 

the  lower-hybrid-drift  instability.  In  order  to  apply  the  results  of  our 

nonlinear  theory  to  the  ionosphere  we  show  in  Fig.  11  a  plot  of  l^i2  “ 

(e*/T. )2(v. /V.. )2(m. /2u  )  versus  n^  ■  9k2  «  9k2p2  after  saturation  for 
l  x  ax  x  e  es 

Vdl/vA  ■  0.56  and  v^ui^  •  0.04.  Since  we  assume  quasi-neutrality  and  the 
ions  are  basically  adiabatic  6n^/n  -  e^/T^),  Fig.  11  can  be  interpreted  as 
an  instantaneous  power  spectrum  of  density  fluctuations  associated  with  the 
lower-hybrid-drift  instability.  The  arrow  at  n2  »  9  denotes  the  fastest 
growing  linear  mode  (for  typical  ionospheric  parameters  B  -  0.3  G,  and 
T^  ~  0.1  ev  this  corresponds  to  X  -  15  cm).  Two  important  aspects  of  this 
figure  bear  comment.  First,  there  is  substantial  power  in  the  long  wave¬ 
length  regime  (n2  <  9)  and  the  spectrum  is  relatively  flat,  consistent  with 
observational  results.  And  second,  the  short  wavelength  spectrum  (n2  >9) 

corresponds  to  a  k  ®  dependence,  somewhat  steeper  than  the  transitional 

—4.8 

regime  which  has  a  k  *  dependence.  Thus,  we  predict  that  the  very  short 
wavelength  power  spectrum  (X  <  15  cm)  has  a  k"6  behavior.  Experimentalists 
are  unable  to  resolve  the  power  spectrum  of  these  very  short  wavelength 
modes  so  that  our  results  cannot  be  verified  at  this  time  but  await  further 
improvements  in  experimental  techniques. 


This  research  has  been  supported  by  0NR  and  NASA.  One  of  us  (J.F.D.) 
was  partially  supported  by  DOE. 


n2  (*lc2  X2) 

The  instenteneous  power  spectrua  is  shown  after  saturation  for 
Vj./v.  ■  0.56  (fn  ■  1.0)  and  ■  0.04  (vn  ■  0.05). 


References 


1*  W.D.  Davis,  A.W.  DeSilva,  W.  Dove,  H.R.  Griem,  N.A.  Krall  and  P.C. 
Liewer,  in  Plasma  Physics  and  Controlled  Nuclear  Fusion  Research 
(International  Atomic  Energy  Agency,  Vienna,  1971),  Vol.  Ill,  p.  289. 

2.  R.C.  Davidson,  N.T.  Gladd  and  Y.  Goren,  Phys.  Fluids  J21_,  992  (1978). 

3.  M.  Tuszewski  and  R.R.  Linford,  Phys.  Fluids  25,  765  (1982). 

4.  J.D.  Huba,  N.T.  Gladd  and  K.  Papadopoulos ,  J.  Geophys.  Res.,  83,  5217 
(1978). 

5.  H.U.  Fahrbach,  W.  Koppendorfer,  M.  Munich,  J.  Neuhauser,  H.  Rohr,  G. 
Schramm,  J.  Sommer  and  E.  Holshauer,  Nucl.  Fusion,  21,  257  (1981). 

6.  R.L.  McPherron,  Rev.  Geophys.  Space  Phys.,  17,  657  (1979). 

7.  N.A.  Krall  and  P.C.  Liewer,  Phys.  Rev.  A,  4 j  2094  (1971). 

8.  R.C.  Davidson,  N.T.  Gladd,  C.S.  Wu  and  J.D.  Huba,  Phys.  Fluids  20,  301 
(1977). 

9.  J.F.  Drake,  J.D.  Huba  and  N.T.  Gladd,  Phys.  Fluids  26,  2247  (1983). 

10.  D.  Winske  and  P.C.  LLewer,  Phys.  Fluids  21,  1017  (1978). 

11.  J.F.  Drake  and  T.T.  Lee,  Phys.  Fluids  24,  1115  (1981). 

12.  R.Z.  Sagdeev  and  A.A.  Galeev,  Nonlinear  Plasma  Theory  (VJ.A.  Benjamin, 
Inc.,  New  York,  1969),  Ch.  II,  p.  37.. 

13.  J.U.  Brackbill,  D.W.  Forslund,  K.  Quest  and  D.  Winske  (private 
communication) . 

14.  R.C.  Davidson,  Phys.  Fluids  21,  1017  (1978). 

15.  J.F.  Drake,  N.T.  Gladd  and  J.D.  Huba,  Phys.  Fluids  24,  78  (1981). 

16.  S.P.  Gary,  Phys.  Fluids  23,  1193  (1980). 

17.  T.H.  Dupree  and  D.J.  Tetreault,  Phys.  Fluids  2J^,  425  (1978). 

18.  P.H.  Diamond  and  J.R.  Myra,  Phys.  Fluids  26,  1481  (1983). 

19.  J.D.  Huba  and  K.  Papadopoulos,  Phys.  Fluids  21,  121  (1978). 

20.  J.F.  Drake,  P.N.  Gusdar  and  J.D.  Huba,  Phys.  Fluids  26,  601  (1983). 

21.  J.D.  Huba  and  S.L.  Ossakow,  J.  Geophys.  Res.,  86,  829  (1981). 

22.  A.  Hasegawa  and  K.  Mima,  Phys.  Fluids  21,  87  (1978). 

23.  W.  Horton,  Phys.  Rev.  Lett.  37^  1269  (1976). 

24.  R.E.  Walts,  Phys.  Fluids  26,  169  (1983). 

25.  P.W.  Terry  and  W.  Horton,  Phys.  Fluids  26,  106  (1983). 


86 


i  -  j  > j;  yy.  • 


H.  Hasegawa  and  D.G.  Pettifor,  Phys.  Rev.  Lett.  50,  130  (1983). 

A.  Hasegawa,  C.G.  Maclennan  and  Y.  Kodarna,  Riys.  Fluids  22,  2122 
(1979). 

D.  Fyfe  and  D.  Montgomery,  Phya.  Fluids  22,  246  (1974). 

D.  Fyfe,  G.  Joyce  and  0.  Montgomery,  J.  Plasma  Phys.  17,  317  (1977). 
S.A.  Orszag,  Stud.  Appl.  Math.  50,  293  (1971). 

H.  Ching,  Phys.  Fluids  16,  130  (1973);  R.E.  Aamodt,  Phys.  Fluids  20, 
960  (1977). 

D.  Winake,  Phys.  Fluids  24,  1069  (1981). 

M.  Tanaka  and  T.  Sato,  J.  Geophys.  Res.  86,  5541  (1981);  Phys.  Rev. 
Lett.  47,  714  (1981). 

Y.-J.  Chen  and  C.K.  Birdsall,  Phys.  Fluids  26,  180  (1983). 

J.D.  Hubs,  J.F.  Drake  and  N.T.  Gladd,  Phys.  Fluids  23,  552  (1980); 
J.R.  Myra,  P.J.  Catto  and  R.E.  Aamondt,  Phys.  Fluids  24,  651  (1981). 
J.D.  Hubs,  P.K.  Chaturvedi,  S.L.  Ossakow,  and  D.M.  Towle,  Geophys. 
Res.  Lett.  5,  695  (1978). 

W.R.  Hoegy,  S.A.  Curtis,  L.H.  Brace,  N.C.  Maynard,  and  R.A.  Heelis, 
Geophys.  Res.  Lett.  9,  993  (1982). 

M.  Singh  and  E.P.  Szuszczewicz ,  to  be  published  in  J.  Geophys.  Res. 
(1983). 

M.C.  Kelley,  R.  Pfaff,  K.D.  Baker,  J.C.  Ulwick,  R.  Livingston,  C. 
Rlno,  and  R.  Tsunoda,  J.  Geophys.  Res.  87,  1575  (1982). 


One  caveat  needs  to  be  added  to  the  discussion  of  the  node  coupling  in 
lower-hybrid-drift  turbulence  presented  in  Sec.  II.  The  dispersion 

A  A 

relation  in  Eq.  (17)  was  derived  by  keeping  only  the  interaction  of  (a>,  k) 

.  -  - 

with  its  nearest  neighbors  (w+,  kj_).  Unfortunately,  the  coupling  to  higher 
order  modes  k  +  pk^  with  p  >  I  cannot  be  neglected  in  general.  This  can  be 
most  easily  understood  in  the  limit  u  »  co+  [see  Eq.  (18)].  Since  we 
have  made  no  a  priori  assumption  on  k,  the  dispersion  relation  should 
remain  invariant  under  the  operation  k  -►  k  +  pk^.  The  dispersion  relation 
in  Eq.  (17)  does  not  have  this  property  and  therefore  cannot  be  correct. 
In  this  limit  there  is  no  small  parameter  which  allows  the  coupling  to 
higher  order  modes  to  be  neglected.  The  coupling  can  be  cut  off,  however, 
by  adding  dissipation  to  the  system.  When  the  modes  k  +  pk^  (with  p  >  1) 
are  heavily  damped,  the  coupling  to  these  modes  can  be  neglected  and  the 
dispersion  relation  in  Eq.  (17)  is  correct.  For  the  parameters  considered 
in  this  paper,  these  modes  are  indeed  heavily  damped  so  that  the  dispersion 
relation  in  (17)  can  be  justified. 
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